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1 Introduction 



It grows the evidence that the universe is undergoing a phase of the accel- 
erated expansion at present epoch. The indications to cosmic acceleration 
appeared not only from the high redshift surveys of type la supernovae [1] 
but also from the anisotropy power spectrum of the Cosmic Microwave Back- 
ground (CMB) [2]. The favored explanation for this behavior is that the uni- 
verse is presently dominated by some form of dark energy. However, neither 
of existing dark energy models is completely satisfactory. Moreover it is very 
hard to construct the theoretical basis for the origin of this exotic matter, 
which is seen precisely at current epoch when one needs the source for the 
cosmic acceleration. 

In the recent papers [3, 4] (see also [5, 6]) it was suggested the gravi- 
tational alternative for the dark energy. It goes as following: the standard 
Einstein action is modified at low curvature by the terms which dominate 
precisely at low curvature. The simplest possibility of this sort is 1/R term, 
other negative powers of the curvature may be introduced as well. Moreover, 
more complicated terms may be suggested like \/{UR + constant) etc. The 
only condition is to dominate over R and to produce the cosmic acceleration 
consistent with current astrophysical data. The interesting feature of the the- 
ories with negative powers of curvature is that they may be expected from 
some time-dependent compactification of string/M-theory as it was demon- 
strated in [7]. Moreover, quantum fluctuations in nearly flat spacetime may 
induce such terms in the same way as the expansion of the effective action 
at large curvature predicts the terms with positive power of (non-local com- 
binations of) curvature invariants (for a recent review, see [8]). 

Clearly, having the gravitational foundation for description of current 
cosmic acceleration seems to be much more natural than the introduction 
by hands of the mysterious dark energy, cosmic fluid with negative pres- 
sure. However, as any other theory pretending to describe current FRW 
universe such a modified gravitational theory should pass number of consis- 
tency checks to be considered as the realistic theory. For instance, as any 
other higher derivative theory, the theory with 1/R may develop the insta- 
bilities [9] (see, however, [10] and discussion in sixth section). As the cosmic 
acceleration from such theory is also instable [3] the question appears: which 
of the instabilities is more realistic or some more modifications of the the- 
ory may be required? From another point, the gravitational action which is 
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the function of only curvature may be presented in a number of ways in the 
equivalent form as scalar-tensor theory with one (or several) scalar (s). Then 
fundamental question is: which action corresponds to our physical world? 
For instance, it has been mentioned [11] that for the model [3] the equivalent 
BD action is not acceptable as it is ruled out by the solar system experiments. 
This indicates against the consideration of the initial modified gravity with 
1/R term as physical theory. In such circumstances it may be a good idea to 
try to search for other variants of modified gravitational theory which still 
excludes the need for the dark energy. Moreover, one can try to construct 
the theory which predicts the inflation at very early Universe and cosmic 
acceleration currently in the combined set-up. 

In the present paper we suggest the new model of modified gravity which 
contains positive and negatives powers of the curvature. Symbolically, the 
lagrangian looks like C = R+ R m + 1/R n where n, m are positive (not neces- 
sary integer ) numbers. The theory may be presented in the equivalent form 
as some scalar-tensor gravity. At large curvature, the terms of the sort R m 
dominate. If 1 < m < 2 the power law inflation occurs at early times. If 
m = 2 the anomaly driven (Starobinsky) inflation occurs at early times. At 
intermediate region, the theory is Einstein gravity. Currently, at low curva- 
ture the terms of the sort 1/R n dominate. These terms serve as gravitational 
alternative to dark energy and produce the cosmic acceleration. It is remark- 
able that such modified gravity does not suffer from the instabilities pointed 
out for the version with the lagrangian L = R + 1/R. Moreover, it may pass 
the solar system tests for scalar-tensor gravity. 

The paper is organized as follows. In the next section we discuss var- 
ious forms of the action for modified gravity, both in Jordan and in Ein- 
stein frames. It is shown also that such theory may have the origin in the 
braneworld scenario by the fine-tuning of surface counterterms. Section three 
is devoted to the study of simplest deSitter solutions for modified gravity. 
The occurrence of two deSitter phases is mentioned. The properties of the 
scalar potential in the scalar-tensor formulation of theory where scalar should 
be identified with the curvature are investigated. In section four, FRW cos- 
mology for such a model (with or without matter) is discussed. It is shown 
that the model naturally admits the unification of the inflation at early times 
and of the cosmic acceleration at late times. In section six we demonstrate 
that higher derivative terms make the dangerous instabilities of original 1/R 
theory to become much less essential at cosmological scales. Moreover, the 
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model easily passes solar system constraints to scalar-tensor gravity. Some 
summary and outlook is presented at final section. 



2 The actions for the modified gravitational 
theory 

Let us start from the rather general 4-dimensional action: 

S = ^J d'x^f(R) . (1) 

Here R is the scalar curvature and f(R) is some arbitrary function. In- 
troducing the auxiliary fields A and B, one may rewrite the action (1) as 
following: 

S = 1 J d 4 xV^ {B {R-A) + f(A)} . (2) 

By the variation over B, A = R follows. Substituting it into (2), the action 
(1) can be reproduced. Making the variation with respect to A first, we 
obtain 

B = f(A) , (3) 
which may be solved with respect to A as 

A = g(B) . (4) 

Eliminating A in (2) with help of (4), we obtain 

S = 1 J d 4 x^ {B (R - g{B)) + / (g(B))} . (5) 

Instead of A, one may eliminate B by using (3) and obtain 

S = ±J d 4 xV^{f'(A) (R-A) + f(A)} . (6) 

At least classically, the two actions (5) and (6) are equivalent to each other. 

The action (5) or (6) may be called the Jordan frame action with auxiliary 
fields. More convenient Einstein frame theory may be worked out as well. 
Under the conformal transformation 

9iiv -> e a g IIU , (7) 
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(i-dimensional scalar curvature is transformed as 



-> e"* (V> - (d - l)Da - {d ~ 1)( f~ 2) g^d,ad v( ^j 



(8) 



Then for d = 4, by choosing 

(7 = -In /'(A), (9) 

the action (6) is rewritten as 



1 f + r- J R 3 
d Xyf—g \R — 



-2 



V/'(A)J 5 P ° f'(A) + f(A) 2 j 



(10) 

l - 2 J d 4 x^9 (R ~ \g p,J d p ad c a - V(a)^j . (11) 



2 

Using a = — In /'(A) = — In 5, the Einstein frame action looks like 
Se — — 

K 

Here 

V(a) = e°g (e"*) - **/ ( 9 (e"")) = ^ - ^ . (12) 

This is the standard form of the scalar-tensor theories where scalar field is 
fictitious one. Although Se = S, we denote the action given in the Einstein 
frame by Se- In (9), it is assumed f'(A) > 0. Even if f'(A) < 0, a = 
— In \f (A) | may be defined. Then the sign in front of the scalar curvature 
becomes negative. In other words, anti-gravity could be generated. 
As a specific choice, we consider 

f(R) = R ~ {R \ iT + b(R - A 2 )" • (13) 

Here we assume the coefficients n,m,a,b > but n, m may be fractional. 
Two last terms may be changed by the sum which includes the terms with 
various negative and positive powers of the curvature: 

f(R)=J2a n R n - (14) 
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Here n can run from negative to positive values. In general, n needs not to 
be necessary integer. It may be also (positive or negative) irrational number, 
that is, F(R) can be rather arbitrary function. 

When first term is absent, n — m > 1, —a = b and A 1 = A 2 the evi- 
dent duality symmetry appears. At very high curvature and supposing both 
cosmological constants to be small, the higher derivative terms dominate. 
Subsequently, with the decrease of the curvature the Einstein gravity domi- 
nates. For low curvature (depending on the choice of Ai) the negative power 
of the curvature may become the leading contribution to the theory. 

If A or R is large 

e" CT = B ~ bm {A - A 2 ) m ' 1 . (15) 
Then if m > 1 and e~ a — > +oo 

V(a) ~ -(bm)~^ (l + -j-^ e^r CT . (16) 

On the other hand when A — Ai or R — Ai is small, we obtain 

e~ a -aniA-A!)- 71 - 1 . (17) 

When e~ a — > +oo 

V(a) ~ -(an)^ (l - e^ a . (18) 

The above arguments indicate that V(a) is not a single- valued function of 
a but at least, there are two branches for m > 2. The above modified 
gravitational action will be our starting point in the attempt to construct 
the Universe where both phases: early time inflation and late time cosmic 
acceleration occur. 



3 The properties of the scalar potential 

Let us discuss the properties of the scalar potential for F(R) given by (13). 
In this section, we concentrate on the case that the matter contribution can 
be neglected. 
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With no matter and for the Ricci tensor R^ u being covariantly constant, 
the equation of motion corresponding to the action (1) is: 

= 2f(R) - Rf'(R) , (19) 

which is the algebraic equation with respect to R. For the action (13) 

= -R+ ( fl *ffin + (m- 2)6 (R - A 2 )™ • (20) 
Especially when n = 1 and m — 2, one gets 

R = R± = WAjTjg (2l) 

If a > 0, one solution corresponds to deSitter space and another to anti- 

A 2 

deSitter. If — < a < and Ai > 0, both of solutions express the deSitter 
space. We may consider other cases, like n = 1 and m — 3. For simplicity, it 
is chosen Ai = A 2 = 0. Then Eq.(20) becomes 

= -R+^ + bR 3 . (22) 
R 

For positive a and b, Eq.(22) has two solutions in general. If we further 
assume a and b is small, there is a solution with large Ri and the one with 
small R s : 

Ri ~ H , i? s ~ v 7 ^ . (23) 

Since the square root of the scalar curvature corresponds to the expansion 
rate of the deSitter universe, the inflation may be generated by a solution 
Ri and the present cosmic speed-up may be generated by R s . One arrives 
at very interesting picture of the universe evolution where modification of 
gravity at high and low curvatures predicts both (early and late time) phases 
of the accelerated expansion. 

We now investigate the potential (12) for the action (13): 

+ b {(m - I) A + A 2 } (A - A 2 ) m ~ 1 
V{A) = - 2 . (24) 



{l+( Iz ^ + 6m(A-A 2 ) m - 1 }' 
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As it is clear from Eq.(2) that A = R, if we regard the scalar curvature R in 
the action (1) with the physical one, A is nothing but the physical curvature 
itself (after equation of motion for A is satisfied) . 
When A — > ±00 

V(A) - ^V m+2 • (25) 

Therefore if m < 2, \/(A) -> 00, if m = 2, \/(A) -> ±, and if m > 0, 
V(A) — > 0. On the other hand, when A — > Ai, one finds 

V(A) — > — (A — Ai) n+1 — > . (26) 
an 

Since the potential (24) is rather complicated, in the following, we only 
consider the case that n = 1, m = 2, Ai = A 2 = 0: 

22 _L_ /,/| 2 

V(A) = , . (27) 

(l + ^ + 26A) 2 

When A -»■ ±00, V(A) -> ^. V(A) vanishes at 

A-0, (28, 



The denominator of the potential (27) vanishes at 

A = A = ~ + a + + a_ , (29) 
06 



± " 2 I 5463 26 26 V V 186 2 A 546 2 
When A — * A , the potential diverges. We should note A < if a, 6 < 0. 

V(-A) has stationary points, where V'(A) = 0, at 

A=(^ , ±v^. (31) 
7 



Then by summarizing the above analysis, we find the rough shape of the 

potential: The asymptotic value of the potential is i. The potential vanishes 

i 

twice when A vanishes and A takes a nee; ative value A = - (2a) 3 in (28). 
The potential is singular when A takes a negative value A = A in (29). The 

potential has three extrema. Two of them are given by positive A: A = (f ) 3 , 
^/3a and one of the three extrema is given by negative A: A = — y/3a. 
When A is small, (9) shows that 

e- ff = ^ + 2^^, (32) 

that is 

A = v^e5 CT . (33) 

Since A — > corresponds to a — > — oo one may only consider the region where 
A and also physical scalar curvature are positive. Then the singularity at 
A = Ao (29) does not appear. 

Since the potential (27) has the following form when A is small, 

V(A) ~ \tf ~ ^=e§ CT , (34) 

in the region where A is small, the universe evolves as the power law infla- 
tion, as we will see later (54). Eq.(34) shows that the potential is increasing 
function for small A. Since the potential is almost flat but monotonically 
increasing for large A, the following cosmological scenario can be consid- 
ered: Let us assume that the universe starts from large but finite positive 
A. Since A corresponds to the physical scalar curvature R, the universe ex- 
pands rapidly because \/R is the expansion rate of the deSitter space. After 
that A rolls down to the small value and the rate of the expansion becomes 
small. When A reaches the small A region, universe begins the power law 
expansion. 

We now give some typical shapes of V(A). In Fig.l, the shape of the 
potential V(A)b versus Ab for ab 2 = ^ is given. When A — > ±oo, V(A) 
approaches to the constant value As mentioned in (28), V(A) vanishes 

at A = and A = - (^)* < 0. When A -> A < in (4), the potential 
diverges. The negative A minimum corresponds to A = —\^3a. In Fig. (2), 
the region A ~ 0.5 is considered. Then two extrema occur, which correspond 
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V(A)b 

5 ■ 
4 ■ 
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Figure 3: V(A)b versus Ab for ab 2 = \. 

to A = (|)* and A = V^a ( (f)* < V^a). If the universe starts at A = 
y/3a, which is locally stable, the universe is in deSitter phase. If by some 
mechanism like the thermal fluctuations etc., A becomes smaller, the universe 

evolves to the power law inflation. The universe may start with A = (f) 3 , 
where the universe is unstable. We should note that the potential with 
ab 2 — ^ is bounded below. In Fig.3, the shape of the potential V(A)b versus 
Ab for ab 2 — | is given. The behavior when A > is qualitatively not so 
changed from ab 2 = case. In the region A < 0, however, the potential 
becomes unbounded below. When ab 2 = the potential is bounded below 
since Aq < —\^3a. Note that A = —\^3a corresponds to the minimum. 
When ab 2 — |, we have A < —\^3a and A = —\/3a corresponds to the 
maximum, then the potential becomes unbounded below. Note once more 
that the above results have physical meaning in terms of our modified gravity 
with negative powers of curvature only when equation of motion for the field 
A itself is satisfied. (Only in this case A may be identified with the curvature 
of the physical metric). 

Thus, we discussed the properties of the scalar potential corresponding 
to quite simple version of the modified gravity. Similarly, the properties of 
other models of such sort may be analyzed. 
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4 FRW cosmology in modified gravity 



Some simple properties of FRW cosmology in modified gravity may be easily 
addressed. One may add the matter to the action (1) with the matter action 
denoted by «S( m ). Then the energy-momentum tensor T( m )^ v can be defined 
by 

2 5S (m) 

T (-W - -^j— ■ (35) 

After the rescaling the metric (7), the energy- momentum tensor T^ m E)^u m 
the Einstein frame (11) is related with T(, ni)iiiU by 

T( m E)nv = e CT T( m ) M1 , . (36) 

Defining the matter energy density p( m ) an d the pressure p( m ) by 

T( m )oo = —P(m)9oo , T( m )ij = P( m )gij , (37) 

the corresponding quantities P( m E) and P( m E) m the Einstein frame (11) are 
given by 

P(mE) = e a p( m ) , P(mE) = e ° V(m) ■ (38) 

We now assume the metric in the physical (Jordan) frame is given in the 
FRW form: 

3 

ds 2 = -dt 2 + d{tf Oijd^dx 3 . (39) 

Here c/ij is the metric of the Einstein manifold, which is defined by the Ricci 
tensor Rij constructed from §ij by Rij = kcjij with a constant k. From the 
conservation of the energy-momentum tensor V^T^ = one gets 

= p(m) + 3if (p( m ) + P(m)) ■ (40) 

Here 

HA, (41) 

and it is assumed p( m ), P( m ) an d also a does not depend on the spatial 
coordinates but on the time coordinate. The matter is chosen to be the 
perfect fluid, which satisfies 

P(m) = WP{m) ■ (42) 
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Eq.(40) can be solved in the usual way 

P{m) = Ca~ 3 ^ , (43) 

with a constant of the integration C. In the Einstein frame, from (38), we 
find 

P{mE) = Ca E - 3 ^e-^ . (44) 
Here cle is the scale factor in the Einstein frame 

cle = e~%a . (45) 

FRW equation in the Einstein frame has the following form: 

3A; k 2 

3#I + = y (PifrES) + P(mE)) ■ (46) 



Here 



H E ^ (47) 
a E 



and p( a E) expresses the contribution from the a field. 

1 /3 



pw = - 2 + ) ■ («) 

With the Einstein frame metric denoted as #£ M „, one obtains 



<rw— • (49) 



V~9e 8a V~9e Sg^ 2 \ g E 

Further using (37), (42), (43), and (45), we find 

1 5S W _ 3W - 1 n ^2. _ 3W- 1 _ 3(1+to) 5^ 

Then the equation of motion for a in the Einstein frame, corresponding to 
the action Se + S'(m) (-Sb is given in (11)) has the following form: 

= 3 (a + 3# £( t) + V'(a) - ^zl^ C d E ' 3{1+w) e-^ . (51) 
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The question now is to solve the system of equations (46) and (51). Before 
going to consider the coupling with matter, we first consider the case of the 
vacuum where C — 0. Especially when the potential is given by (34) as 
A = R is small, a solution is given by 

aE = a E o[—J , cr = --In— . (52) 

Here t E is the time coordinate in the Euclidean frame, which is related with 
the time coordinate t in the (physical) Jordan frame by e^dts = dt. As a 
result 

3*1 = t , (53) 
and even in the physical (Jordan) frame the power law inflation occurs 

a = e^chE oc t E oc t 2 , (54) 

which is consistent with the result in [3]. 3 Hence, at small curvature the 
(instable) cosmic acceleration is predicted by the terms containing inverse 
curvature. If the present universe corresponds to the above power law infla- 
tion, the curvature of the present universe should be small compared with 
that of the deSitter universe solution in (21) with Ai = 0. As the Hub- 
ble constant in the present universe is (lCT 33 eV) \ the parameter a, which 
corresponds to /i 4 in [3], should be much larger than (10 _33 eV) 4 . 

One may also consider the case that the curvature is large. When A or 
the curvature is large, the potential (27) becomes a constant: V(A) — > ^. 
Then from (51), we may assume, if we neglect the contribution from matter, 
that 4> is small. From (46) it follows 

3 ^ + 2^~S' (55) 

which shows that the spacetime is deSitter. For the case of A; = 0, H E 
becomes a constant: 



in [3] can be identified as a — — 
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Then universe expands exponentially. 

The matter contribution to the energy- momentum tensor may be ac- 
counted for in the same way. When it is dominant compared with one from 
cr, Eqs.(46) and (51) can be reduced as 

= SHl + J^-^CdE-^e-** 1 ' , (57) 

= 3(a + 3H E &)-^^K 2 Ca E - 3 ^e- 3 ^ . (58) 
When w — |, which corresponds to the radiation, from (58), one finds 

Here d and a E0 are constants. Eq.(59) expresses the redshift of a. When 
w = 0, which corresponds to the dust, if we can regard a is almost constant, 
we find 

a = a (^j 3 . - ^ 2( *~ t0 W 3(1+w) e-^ (60) 

The obtained results are not changed from those in [3] where the possibility 
of cosmic acceleration in 1/R model was established. 

When n = 1 or in general n is an odd integer, as is clear from (18), A = 
is not a minimum nor local minimum although V'(A) = at A = 0. Of 
course, since A — > corresponds to a — > +oo, it might be unnecessary to 
consider the region A < 0. In the region A > 0, A = can give a minimum 
of the potential. On the other hand, when n is an even integer, A = is 
at least local minimum. Then if inflation occurs in the region A ~ 0, the 
inflation should be stable. 

We now consider the general case that f(A) is given by (13). When R = 
A ~ Ai is small, the potential is given by (18). Neglecting the contribution 
from the matter fields, solving (46) and (51), we obtain, instead of (52), 

(n+l)(2n+l) 

. ft E \ W 2 w+1 ) t E 

aE = aE «\Tj ' a = -^TW ln ^' (61) 

Instead of (53) the physical time is 

(n + 2)tf* = t . (62) 
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Then the power law cosmic acceleration occurs in the physical (Jordan) 
frame: 

A (n + l)(2n + l) 

a oc t ™+ 2 . (63) 

It is quite remarkable that actually any negative power of the curvature 
supports the cosmic acceleration. This gives the freedom in modification of 
the model to achieve the consistency with experimental tests of newtonian 
gravity. 

On the other hand, if the scalar curvature R = A is large, the potential 
is given by (16). Solving (46) and (51) again, one obtains 

(m-l)(2m-l) 

A ft E \ (—2)" 2(m-l), t E , x — K 
a E = aEo(-=-) , <* = -4 4~ m 4 ? (m -2)t E m ~ a = t . 

(64) 

Thus in the physical (Jordan) frame, the universe shrinks with the power law 
if m > 2 4 : 

(m-l)(2m-l) 

a oc t m - 2 . (65) 

Of course, if we change the arrow of time by t — > t — t, the inflation occurs 
with the inverse power law and at t — to, the size of the universe diverges. 
It is remarkable that when m is fractional (or irrational) and 1 < m < 2, 
Eqs.(64) and (65) are still valid. Then the power in Eq.(65) becomes positive, 
the universe evolves with the (fractional) power law expansion. 

Anyway, it is interesting that our model may unify both phases: early 
time inflation (for m < 2 or m = 2) and current cosmic acceleration. 

In [3], from the dimensional analysis it is predicted for /z 4 = a (a for 
n = 1 case) to be \x ~ H ~ 10~ 33 eV. Here H Q is the Hubble constant in the 
present universe. From the dimensional analysis, one may choose 

a= (lO" 33 eV) 2(n+1) , b= (l(T 33 eV) 2( - m+1) . (66) 

In principle, these parameters should not be small ones, unless this is pre- 
dicted by experimental constraints. 

Adding the (dominant) matter contribution to the energy-momentum 
tensor, the situation is not changed from n — 1 and m — 1 case and the case 
of [3]. We obtain (59) and (60) again. Thus, the possibility of cosmic accel- 
eration in the model with negative powers of the curvature is demonstrated. 

4 For Tji = 2 the well-known anomaly driven (Starobinsky) inflation occurs. 
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Moreover, terms with positive powers of the curvature in such a model may 
realize the inflation at early times. 



5 Simplest tests for modified gravity 

In this section, we discuss the (in) stability of our higher derivative model 
under the perturbations. Simplest constraint to the theory parameters from 
the equivalent BD-type gravity is also analyzed. 

In [9], small gravitational object like the Earth or the Sun in the model [3] 
is considered. It has been shown that the system quickly becomes instable. 

The general 5 equation of motion corresponding the action (1) with matter 
is given by 

\g, v f{R) - R^f(R) - g^Of'(R) + V,V v f{R) = -yT (m)/lv . (67) 

Here T( m )^ v is the energy-momentum tensor of the matter. By multiplying 
to (67), one arrives at 6 



K 



2 



nR + fW) wpRWPR + WW) ~ WW) " WW) T ' (68) 

Here T = T^ p . Let f(R) is given by (13). Then in case of the Einstein 
gravity, where a = b = 0, the solution of Eq.(68) is given by 

R = R = ~T . (69) 

The perturbation around the solution (69) may be addressed 

R = R + R 1 , flflxl < \Rq\) . (70) 
Then by linearizing (68), we obtain 

= OR + l^mv p R WR 0+ f '^ 2fW r ° 



fW(R ) p u u 3/W (Ro) 3/W (R ) 3/( 2 )(i? ) 

f (3) (Rn) 

+ Di? i + 2 ri 2> D W^i + U(R )Ri ■ (71) 



5 We consider the case that the Ricci tensor is not covariantly constant 

6 Thc convention of the signature of the spacetime here is different from that in [9]. 
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Here 



f( 1 \R Q )f^(R )R / (1) (^o) , 2/( J R )/ (3) (^o) Rof® 



3/ (2 )(jRo) 2 /(2)(jRo) 3/( 2 )(i? ) 2 /< 2 >(i2o) 



.(72) 



If U (Rq) is negative, the perturbation Ri grows up exponentially with time. 
The system becomes instable. In the following, we neglect the terms of the 
sort Vp-Ro- In case of the model [3], where b = Ai = A 2 = 

U(Ro) = -R + ^. (73) 

In order to describe the universe acceleration in the present epoch, we may 
have [3] 

/i^ 1 = cT3 ~ 10 18 sec ~ (KT 33 eV) _1 . (74) 

(As discussed after Eq.(54), if the present universe corresponds to the power 
law inflation, a can be much larger than (10~ 33 eV) 4 .) By using the estima- 
tions in [9], 



R 3 



(lO^secW-^O , R ~ (I0 3 sec)- 2 ( , (75) 
\gcm~ J / ' ygcm - - 3 / 

we find the second term in (73) dominates. Here p m is the mass density of the 
gravitating body. Since R is negative from (69), U(R ) becomes negative 
and the microscopic gravitational instability occurs. 

First we consider the case that n — 1, m — 2, and Ai = A 2 = 0. For 

b » r^i , (76) 

one gets 

U(Ro) ~ y < • (77) 
Then the instability shows up again. We should note, however, that since 
llM ^ la fr° m (75), the (macroscopic) instability development takes quite 
long time. In fact, the time for instability to occur is significally improved 
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(by the order of 10 29 ). Clearly, other higher derivative terms may significally 
improve the estimation as is shown below. 

Eqs.(74), (75), and (76) indicate that b' 1 < (10 n eV) 2 if we assume 



p m ~ lg cm -3 . Then ^ <C (1CT 17 ) 2 , which does not seem to be consis- 
tent with the bounds obtained from the observation [12] (see also [13] where 
HD inflation was discussed). The action of model in [13, 12] contains the 
Einstein-Hilbert term and the i? 2 -term. There are two ways to overcome 
this bound. First of all, one can include the contributions from other terms 
containing some powers of the curvatures and/or the contribution from the 
matter fields and the quantum effects like conformal anomaly during the in- 
flation. For instance, account of conformal anomaly gives the way for trace 
anomaly driven inflation. The bound might be changed. From another side, 
as it follows from the discussion after Eq.(54), if the present accelerating uni- 
verse corresponds to the power law inflation, a should be much larger than 
(lCT 33 eV) 4 . If we choose a to be large enough in (73) so that a ^> Rq, the 
first term in (73), which is positive, dominates. Then the instability [9] does 
not exist from the very beginning. 

For n — 1, m — 3, and Ai = A 2 = one can assume 

b » > (78) 



and 



Then if 



CD 1 



ft< 5SI- (80) 

U(Ro) is positive and the system becomes stable. From (75), we find 
jp. Then the condition (80) is compatible with the assumption (78). Thus, 
the addition of the terms with the positive powers of the curvature may 
salvage the modified gravity as less instabilities occur or their development 
is macroscopic and other effects may prevent them. 

It has been mentioned in ref. [11] that 1/R model which is equivalent to 
some scalar-tensor gravity is ruled out as realistic theory due to the con- 
straints to such theories. Let us study if our theory may pass these con- 
straints. The coupling of the a-field with matter is always of the same order 
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with gravity [14]. Then if a field has a small mass, the present model cannot 
be realistic. One may consider the case that the present universe corresponds 
to one of the solutions (21) with A 1 = 0, which corresponds to a minimum 
in (31) of the potential V(A): R = A = ^/3a. We now calculate the square 
of scalar mass, which is proportional to V"(a). When A = \^3a 



d 2 V(a) 



da 2 



A=V3a 




d 2 V(A) \ 
dA 2 J 



3a (| + 2bV3a) 



2 



(81) 



Choosing 



the mass of a becomes large, 
system tests [17]. 



(82) 

Thus, HD term may help to pass the solar 



6 Discussion 

In summary, we considered the modified gravity which naturally unifies two 
expansion phases of the universe: inflation at early times and cosmic acceler- 
ation at current epoch. The higher derivatives terms with positive power of 
the curvature are dominant at the early universe providing the inflationary 
stage. The terms with negative power of the curvature serve as gravitational 
alternative for the dark energy making possible the cosmic speed-up. It is 
shown that such a theory is stable unlike to its more simple cousin with only 
Einstein and 1/R terms. Moreover, it may pass the simplest solar system 
constraint for scalar-tensor gravity equivalent to original modified gravity. 
Clearly, more checks of this theory should be done in order to conclude if 
the model is realistic one or not. In any case, there is some space for modi- 
fications by adding the terms with other powers of the curvature. Moreover, 
the use of other curvature invariants (Ricci and Riemann tensors) may be 
considered as well. 

It seems very attractive that modified gravity which may origin from 
string/M-theory eliminates the need for dark energy. Clearly, deeper inves- 
tigation in this extremely interesting direction connecting the cosmological 
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constant problem, early time inflationary and current accelerating universe is 
necessary. It is a challenge to find the solution for above fundamental cosmo- 
logical problems from the first principles. The search for true gravitational 
action may be the step in this direction. 
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A Modification of the gravity from the infla- 
tionary potentials 

It is very interesting that modified gravity may be often predicted by the 
standard inflationary cosmology with inflaton field if the corresponding po- 
tential is not trivial. Let us consider the model coupled with the scalar field 
with exponential potential: 

S = J d'x^g i-^-R - -sTd^t - V e^\ . (83) 

Here M P is the 4-dimensional Planck mass. It has been shown [15] that if 
a < \f2, the potential becomes shallow and supports the inflation. Recently 
there appeared some arguments towards to possibility of the inflation and 
cosmic acceleration in the theory based on the braneworld scenario with 
similar exponential potentials [16]. If the metric is rescaled 

9,u - ei^ v , (84) 
the action (83) can be rewritten as 

S = f d^x^-g (¥k e *yfi+R - V eM"+Vl)^ . (85) 
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Since the kinetic term for vanishes, we may regard to be an auxiliary 
field. By the variation over 0, one gets 







M P 4/5 



4V37T 



Vi*R-3-Vto( a + */T\ eM^Vl) 4- P > (86) 
Mp \ V 3 I 



which may be solved as 



STrv/6 (a + Vo 
By substituting (87) into the action (85), we obtain 



R 



d = r- 1 — / a XyJ—g 



MlR 



8nV6 ( 



a 



(87) 



(88) 



Then the action (83) is equivalent to the modified gravity with the fractional 

power °f ^ ne scalar curvature. If — < a < — , the power of the 

curvature becomes negative again. 

Similarly, one can consider other inflationary potentials. It is clear that 
many of them lead to modification of the gravitational theory , supporting 
the point of view that inflaton may be not physical but rather some fictitious 
scalar field. 
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